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Abstract. We represent the generators of the SU(iV) algebra as bilinear combinations of Fermi operators 
with imaginary chemical potential. The distribution function, consisting of a minimal set of discrete imag- 
inary chemical potentials, is introduced to satisfy the local constraints. This representation leads to the 
conventional temperature diagram technique with standard Feynman codex, except that the Matsubara 
frequencies are determined by neither integer nor half-integer numbers. The real-time Schwinger-Keldysh 
formalism is formulated in the framework of complex equilibrium distribution functions for auxiliary semi- 
fermionic fields. We discuss the continuous large N and SU(2) large spin limits. We illustrate the application 
of this technique for magnetic and spin-liquid states of the Heisenberg model. 

PACS. 7 5.10.-b, 75.10.Jm, 75.40.Gb, 71.10.Fd 



Several approaches have been proposed for the descrip- 
tion of spin systems in statistical physics. Methods of func- 
tional integration based on various representation of spin 
operators such as Fermi, Bosc, Majorana, supersymmetric 
or Hubbard operators have been 

applied to man y p roblems involving quantum spins and 
pseudospins |L2 



13 



y,y,y,y,|ij,y,grjy,^,^3|. The dif- 

Acuity with the representation of spin operators is con- 
nected with the fact that spins possess neither Fermi nor 
Bosc statistics. The commutation relations for spins are 
determined by the SU(2) algebra, leading to the absence 
of a Wick theorem for SU(2) generators. The Gaudin pi| 
theorem existing instead makes it impossible to construct 
a simple diagram technique directly for spin operators. To 



resolve this problem, various representations 
^,P, ^0| , [lll] have been introduced. Nevertheless 



HUM 

the repre- 



X 



sentation of spins as a bilinear combination of Fermi/Bose 
operators enlarges the dimensionality of Hilbert space where 
these operators act. Thus, the spurious (unphysical) states 
should be excluded from the consideration resulting in 
a constraint requirement. Basically, different representa- 
tions cure the constraint problem in a different way. Nev- 
ertheless, the usual price for simplicity is the replacement 
of the local constraint on each point containing the spins 
by a so-called global constraint, so that the restriction is 
fulfilled only in the average over all sites, ft is known that 
such a replacement results in uncontrollable approxima- 
tions for quantum spins (especially in low dimensions). 
Although the use of a global constraint is questionable 
for SU(2) systems, it becomes more reasonable for higher 
SU(AT) groups, especially in the "large N limit". The cor- 
responding approach is known as "1/N expansion", jl7], 
[l8| , 19 ,[20 , 2lJ successfully describing the strong coupling 



limits of the Kondo imp urity |i2| , Anderson lattice [[13| 
[li], [[5) and Hubbard |lj| models and also SU(AQ Heiscn- 
berg antiferromagnets on a square lattice 00^1^113, 
p2[ shedding light on the mechanism of high T c supercon- 
ductivity in cuprate compounds. Although SU(iV = 2) 
models are of primary physical interest, the SU(iV ^ 2) 
models can be considered as " approximate models" where 
an "exact solution" can be gained in contrast to "exact 
models" where the "approximate solution" is hard to ob- 
tain |^,|l). The simplification arises due to expansion in 
the inverse number of "flavors" 1/N, making it possible 
to start with mean-field solution and systematically find 
corrections to it. 

The goal of this paper is to consider a semi-fermionic 
representation for SU(JV) generators for arbitrary (not 
necessary large) N, applying a different idea of constraint 
realization. This idea is know as Popov- Fedotov |25| rep- 
resentation being initially proposed for S = 1/2 and 5 = 1 
spins. Based on an exact representation of spin operators 
as fermions with imaginary chemical potential, this rep- 
resentation resulted in the conventional Feynman temper- 
ature diagram technique, nevertheless providing a rigor- 
ous treatment of the local constraint. pfj|, p7||2^ ,p9|, |3Q| , |3l|l 
In this paper we give a generalization of this method to 
SU(iV) and we also construct the real-time formalism for 
the semi-fermionic spin representation. 

The SU(iV) algebra is determined by the generators 
obeying the following commutational relations: 



(1) 



where a, j3 = 1 , . . . , N. We adopt the definition of the Car- 
tan algebra §§] of the SU(N) group {H a } = similar 
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n c 



m 



a) b) c) 

Fig. 1. Rectangular Young tableau to denote a SU(iV) repre- 
sentation (a), single column tableau corresponding to n c = 1 
(b) and single row tableau standing for spin S — n c /2 repre- 
sentation of SU(2) group (c). 

to the one used in p7| , noting that the diagonal genera- 
tors S" are not traceless. To ensure a vanishing trace, the 
diagonal generators should only appear in combinations 

JV N 

J2 s *sz ^th j2 Sa = ( 2 ) 

a— 1 a— 1 

which effectively reduces the number of independent diag- 
onal generators to TV — 1 and the total number of SU(7V) 
generators to iV 2 — 1. 

For SU(2) one recognizes the usual spin operators 

5 2 = 5+, Sl = S~, Sl-S\ = 2S z (3) 

with the usual commutation relations |53j and the Pauli 
matrices as generators of the SU(2) group. We shall not 
confine ourself to some special type of Hamiltonian. Never- 
theless, it's worthwhile to mention that the SU(7V) gener- 
alization of the Hciscnberg model is given by the following 
expression |]],|o|,[2l[] 

*=^££W?c*) (4) 

<ij> a(3 

On each site, there may exist many particles, whose sym- 
metry properties define a specific representation of SU(-/V). 
The most transparent way to visualize an irreducible SU(iV) 
representation are Young tableaux (34[|35| . Instead of the 
general Young tableau, specified by N — 1 integers, for ex- 
ample the lengths of the rows, we restrict us for the main 
part of the paper to rectangular tableaux, with 1 < m < 
N rows and n c > 1 columns, illustrated in Figjl]. 

The familiar spin is given by N = 2, m = 1, n c = 25*, 
so the Young tableau contains one row of 25* length, with 
only one box for 5 = 1/2. Also of special importance are 
the tableaux with n c = 1, giving the N — 1 fundamental 
representations of SU(iV). 

The Sp generator may be written as biquadratic form 
in terms of Schwinger boson operators (2^| . |53| : 

S$ = biy* (5) 



and a constraint as follows 

N 

J2 b U baq = s l n c (6) 

a=l 

where p = 1, ...,m is the number of "colors". 

The equivalent fermionic representation of the gener- 
ators of SU(7V)@ is given by 

^ = E c -^ a (7) 

a 

where the "color" index a,b = 1, n c and the n c (n c +l)/2 
constraints 

N 

E c -c ab = ^m (8) 

restrict the Hilbert space to the states with m*n c particles 
and ensure the characteristic symmetry in the color index 
a. The antisymmetric behavior with respect to a is a direct 
consequence of the fermionic representation. 

Let us consider the partition function for the Hamil- 
tonian, expressed in terms of SU(iV) generators 

Z s =Trexp(-0Hs) (9) 

For SU(2) 5 = 1/2 and S = 1 it is possible to map the 
spin partition function onto a fermionic partition function 
where the chemical potential of fermions is purely imagi- 
nary 

Z s = ATr exp {-/3{H F - fiN F )) = AZ F (10) 

with fi = -ittT/2 and A = i nM ° for S = 1/2 and fi = 
—iirT/3 and A — (i/v3)™ s " e for S = 1, n S ite denotes 
the number of sites in a lattice. This results in usual 
Feynman-like diagram technique built up with the help 
of auxiliary Fermi (Grassmann) fields. The corresponding 
Matsubara frequencies for Popov-Fedotov (PF) fermions 
after applying the generalized Grassmann boundary con- 
ditions J25| read as w n = 2irT(n + 1/4) for S = 1/2 and 
u> n = 2irT{n+l/i) for 5=1. The imaginary chemical po- 
tentials are important for the realization of an exclusion 
principle providing the fulfillment of the general identity 

Z s = 7>cxp(-/lff s ) = TreM-0H F )5 nFA (11) 

The Popov-Fedotov representation has been generalized 
for arbitrary values of spin 5 for the SU(2) group in [ |36| by 
introducing the distribution of discrete chemical potentials 
with j being the site index, for PF fermions: 

z s = \{j dn(j)P(nU))Z F (nU)) (12) 

3 

For the SU(iV) algebra we shall try to find the partition 
function in a similar manner 

Zs= [ l[dn(j)P{n(j))Trexp (-p(H F - v(j)n F )) = 
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3 



(13) 



We use the path integral representation of the partition 
function 

Zs/Z° s = jY[d^j)P^(j))exp(A)/ 



l[dfx(j)P{p(j))exp(A ) 



(14) 



where the action A and Ao are determined by 

A = Ao- / (1tH f (t), 
Jo 

A = EE dTa k (j,T)(d T + fi(j))a k (j,r) (15) 
j fe=i Jo 

and the fermionic representation of SU(iV) generators (0) 
is applied. 

To begin with we confine ourselves to two particular 
cases of SU(7V) with n c = 1 (corresponding to an effective 
"spin size" S = 1/2 and in the language of Young tableaux 
described by one column) and SU(2) for arbitrary value 
of n c = 2S (one row Young tableau). 

Let us first consider n c = 1. We denote the corre- 
sponding distribution by PN,m(p(j)), where m is the num- 
ber of particles in the SU(iV) orbital, or in other words, 
1 < m < N labels the different fundamental representa- 
tions of SU(iV). 



N 



(16) 



fc=i 



To satisfy this requirement, the minimal set of chemical 
potentials and the corresponding form of -P/v,m(Mj)) are 
to be derived. 

Let us classify the states in Fock and spin spaces. We 
note that the dimension of the Fock space is dimHp — 2 N 
and spurious states should be excluded. Thus, there are 
v(N, m) = C% 1 = N\/(m\(N — to)!) physical states which 
can be obtained from the vacuum state = |0, ...,0) as 



follows 



N 



<P 



M 

phys 



(IW)*o 



(17) 



i=i 



or from the| 1, 1, 0, 0) state by transferring the oc- 



N — 7 



cupied states from left to the right side using the group 
generators. 

To derive the distribution function we use the following 
identity for constraint ( |l6| ) expressed in terms of Grass- 
mann variables 



<5„ 3 ,m = sin (Tr(n 3 - to)) / sin I — I 



(18) 



Substituting this identity into (|T] 
T4) on gets 



and comparing with 



PN.,m{Kj)) 



where 



1 



JV 
k=l 



exp 



N 



(2k - 1) 8(n(j) - Mfe) 



iirT 



{2k -I) 



(19) 
(20) 



Since the Hamiltonian is symmetric under exchange of 
particles and holes if the sign of the chemical potential is 
changed simultaneously, we can simplify ( |Tg| ) to 



2i W2i 

PN, m (v(j)) = sin 

fc=l 



2k- 1 



(21) 

where [N/2\ denotes the integer part of N/2. As the dis- 
cussion below will show, this is the minimal represen- 
tation of the distribution function corresponding to the 
minimal set of the discrete imaginary chemical potentials. 
Another distributions function different from (|2l]) can be 
constructed when the sum is taken from k = N/2 + 1 to 
N. Nevertheless, this DF is different from (|2l]) only by the 
sign of imaginary chemical potentials ji k = /u£ = — [i k and 
thus is supplementary to (pl|). 

Particularly interesting for even N is the case when the 
SU(iV) orbital is half-filled, to = N/2. Then all chemical 
potentials are weighted with equal weight 



Pn,n/2(p(j)) = I EM)^ 1 ^') - /**) 



N/2 



(22) 



fe=i 



Taking the limit N — > 00 one may replace the summation 
in expression (22) in a suitable way by integration. Note, 
that taking N — > 00 and to — > 00 we nevertheless keep 
the ratio m/N =1/2 fixed. Then, the following limiting 
distribution function can be obtained: 



Pn,n/2(Kj)) 



N- 



JL 

2ni 



exp 



N 



(23) 



resulting in the usual continuous representation of the lo- 
cal constraint for the simplest case n c = 1 (compare it 
with Q)) 

N 

Z s = Tr(exp {-(3H F ) Sfa - -)) (24) 



We note the obvious similarity of the limiting DF (|2£ 
with the Gibbs canonical distribution provided that the 
Wick rotation from the imaginary axis of the chemical 
potential fi to the real axis of energies E is performed and 
thus n(j)N/2 has a meaning of energy. 

Up to now the representation we discussed was purely 
fermionic and expressed in terms of usual Grassmann vari- 
ables when the path integral formalism is applied. The 
only difference from slave fermionic approach is that imag- 
inary chemical potentials are introduced to fulfill the con- 
straint. Nevertheless, by making the replacement 



ak{j,T)) -> a k (j,T) exp 



2k- r 



13 N 
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a>k(j,T) -> a fe (j,T)exp 



iiTT 2k — 1 



(25) 



we are coming to generalized Grassmann (semi-fermionic) 
boundary conditions 



2k -1 



ak(j,(3) = a k (j,0) exp in 



a>k{j, P) = ak(j, 0) exp -in 



N 
2k- 1 



N 



(26) 



This leads to a temperature diagram technique for Green 
functions 



G a0 (j,T)=-(T T a a (j,T)a p (j,O)) 



(27) 



toe** , 




SU(2) 


=2 

-\* 1 








J 




of semi-fermions with Matsubara frequencies different from 

both Fermi and Bose representations. Fi g . 3. Graphical representation of exclusion principle for 

The minimal set of Matsubara frequencies w„/ (2nT) SU(JVj semi-fermionic representation with even N,n c = l (we 

forms for SU(iV) with even N the triangle table shown in use p = inT/2 for SU(2) and m = inT/4, fi 2 = 3inT/4 for 

Figf SU(4)). 



Fermi 




' 2N' 



+ 



2tiT 



n n+1/2 

Fig. 2. The minimal set of Matsubara frequencies for SU(N) 
representation with even N. 



The exclusion principle for this case is illustrated on 
Fig.||, where the first two groups SU(2) and SU(4) are 
shown. The first point to observe is that the spin Hamil- 
tonian does not distinguish the n particle and the n hole 
(or N—n particle) subspace. Due to Eq. ( |20| ) the two phase 
factors exp(/3^n) and exp(/3 fi(N — n)) accompanying these 
subspaces in Eq. ( pl| ) add up to a purely imaginary value 
within the same chemical potential, and the empty and 
the fully occupied states are always canceled. In the case 
of N > 4, where we have multiple chemical potentials, the 
distribution function P(/Lt) linearly combines these imagi- 
nary prefactors to select out the desired physical subspace 
with particle number n = m. 

In Fig.||, we note that on each picture the empty and 
fully occupied states are canceled in their own unit cir- 
cle. For SU(2) there is a unique chemical potential \i = 
±iirT/2 which results in the survival of single occupied 
states. For SU(4) there are two chemical potentials (see 
also Fig.||). The cancellation of single and triple occupied 
states is achieved with the help of proper weights for these 



states in the distribution function whereas the states with 
the occupation number 2 are doubled according to the ex- 
pression (^2]). In general, for SU(iV) group with n c — 1 
there exists N/2 circles providing the realization of the 
exclusion principle. 

We consider now the generalization of the SU(2) alge- 
bra for the case of a large moment S with 25+1 projec- 
tions. Here, the most convenient fermionic representation 
is constructed with the help of a 25 + 1 component Fermi 
field Ofe(j') provided that the generators of SU(2) satisfy 
the following equations 

S-i 

5+ = VS(S + l)-k(k + l)al +1 (j)a k (j) 

k=-S 



S~= E ^(S+1)-*(*-1)4-i(J>*(j) 

k=-S+l 



s z = E fc 4(iK(i) 



(28) 



k=-S 



such that dimHp = 2 2S+1 whereas the constraint reads 
as follows 



k=S 



(29) 



k=-S 



We consider the distribution function for arbitrary I for 
the sake of generality. It describes the orbital part of an 
atomic subshell with orbital quantum number 5, with / 
particles present. We denote the corresponding distribu- 
tion function of the chemical potential by P2s+i,i(fJ-(j))- 
Following the same routine as for SU(iV) generators and 
using the occupancy condition to have I (or 25 + 1 — I) 
states from the 25+1 states filled, one gets the following 
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distribution function, after using the particle-hole sym- 
metry of H$: 



2S+1A 



2i 



25+1 



LS+1/2J 



fc=l 



tt/- 



,2k- 1 
'25+1 



^(M(i)-Mfe) 



(30) 

where the chemical potentials are /x^ = —inT(2k—l)/ (25+ 
1) and fc = 1, [5 + 1/2J, similarly to Eq.@. 

In the particular case of the SU(2) model with I = 1 
for some chosen values of spin 5 the distribution functions 
are determined by the following expressions 



P2,i (m(j)) 



for 5 = 1/2 



7 7 7T7 1 

p 3 ,i(m(j)) - ^, 2 (m(j)) = ^= W) + — ) 



(31) 



(32) 



for 5 = 1. 

This result corresponds to the original Popov-Fedotov 
description restricted to the 5 = 1/2 and 5=1 cases. We 
present as an example some other distribution functions 
obtained according to general scheme considered above: 

p 4 ,iO) =P4, 3 (m) - 

= ^(^+^)+%+^) (33) 
for 5 = 3/2, SU(2) and 

PM = \ (S(M + Sip + ^)) (34) 

for effective spin "5 = 1/2", SU(4), 




1 r/ * 7rr 
7!^ + — 




1 p/ 3i7rT x 



(35) 

for 5 = 2, 5C/(2) etc. 

A limiting distribution function corresponding to Eq. 
( p3| ) for the constraint condition with arbitrary I is given 

by 



S— *oc 



27ri 



exp(-/3//i(j)) 



(36) 



For the case I = m = N/2 -> oo and 5 = (AT - l)/2 -»• oo 
the expression for the limiting DF Poo,i(n(j)) coincides 
with (23). We note that in 5 — > oo (or N — > oo) limit con- 
tinuum chemical potentials play role of additional U(l) 
fluctuating field whereas for finite 5 and N they are char- 
acterized by fixed and discrete values. 

When 5 assumes integer values, the minimal funda- 
mental set of Matsubara frequencies is given by the table 
in Fig.g 

The exclusion principle for SU(2) in the large spin limit 
can be also understood with the help of Fig.|j and Fig.||. 




n + 1 n + | 

, 3 \ Fermi 



n li u 1 n 



2jiT 



1 h 

n n+1/2 

Fig. 4. The minimal set of Matsubara frequencies for SU(2) 
representation for integer values of the spin and 1 = 1. 




Fig. 5. Graphical representation of exclusion principle for 
SU(2) semi-fermionic representation for S = 1 and S — 2. 
For any arbitrary integer value of spin there exists S circle 
diagrams corresponding to the S different chemical potentials 
and providing the realization of the exclusion principle. 



One can see that empty and fully occupied states are can- 
celed in each given circle similarly to even-N SU(iV) al- 
gebra. The particle-hole (PH) symmetry of the represen- 
tation results in an equivalence of single occupied and 25 
occupied states whereas all the other states are canceled 
due to proper weights in the distribution function (|o|). In 
accordance with PH symmetry being preserved for each 
value of the chemical potential all circle diagrams (see 
Fig. 3, Fig. 5) are invariant with respect to simultaneous 
change //<->-// and n partic i e <-> n ho i es . 

Let us make few comments concerning the general rect- 
angular Young tableau of size n c * m. The fermionic repre- 
sentation (fjj) is characterized by an N*n c component field 
with n c identical diagonal constraints and n c (n c — l)/2 
offdiagonal constraints (||). The effective "filling" deter- 
mining the number of fermions on each site is mn c . How- 
ever, not all of these (n c N)\/((n c m)l((N — m)n c )\) states 
are representing proper physical states. One should take 
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into account the constraints Eq. (g) to obtain the complete 
set. The number of physical states of a rectangular Young 
tableau is given by the expression: 



N-m 



v{N, m, n c ) = 



°AT+Ti c -l Ar+n c -2--- L/ m+n c 



N - 



oo 



c. 



t-i (3 



D + OO 





-l°7V+n c -2- 






sym 
-l°m+n c -2- 





(37) 



While the diagonal part of Eq. (Q) could be satisfied with 
color-dependent chemical potentials [i a (j)i coupling only 
to ^2 a cT a c Qa , an exclusion procedure for the off-diagonal 
constraints needs either projection operators or an a priori 
restriction on the trace, using e.g. coherent states Ej. 

Another generalization is applicable for a broader range 
of cases. The general Young tableau (not necessarily rect- 
angular), representing any irreducible representation (p), 
can be described in the context of our approach in the fol- 
lowing way. The generators Sm are expressed as matrices 

(S£ ) )p a = (^ ) \T£^\^) (38) 

with Tm^^ — T^n + Tm being the generators in a suit- 
able direct product of representations (k) and (I) and the 

states |^i p ^) are obtained in terms of Clebsch-Gordon co- 
efficients 

l^ p) > = E( fc > l > * Hp; ")l^ fc) ) x (39) 

This procedure can be iterated until (fc) and (I) are fun- 
damental irreducible representations of SU(iV). The size 

of the matrices S$ is equal to the dimension of the repre- 
sentation, v^. The trace is now easily expressed in terms 
of fermionic fields, enforcing the constraint 8 n .\ with 
the distribution of chemical potentials (see Eq.(]2lj)) 



P v iv) And)) = 



2i 



k=l 



2k -1 



%0') - Wfe) 
(40) 

For the simple case of SU(2), which yields only single-row 
tableaux, this procedure gives the fermionic representa- 
tion described in Eqs.( p8[]30| ). In the case of single-column 
tableaux for SU(iV), however, and in the general case of 
mixed symmetry, it does not fully use the fermionic com- 
mutation properties. Therefore, it is in general not the 
representation with the minimal number of fermions and 
the minimal number of chemical potentials in -P(/i). 

We discuss finally the real-time formalism based on 
the semi-fermionic representation of SU(iV) generators. 
This approach is necessary for treating the systems be- 
ing out of equilibrium, especially for many component 



Fig. 6. The Keldysh contour going from — oo — » oo — > — oo 
in real time. The boundary conditions on the imaginary time 
segment determine the generalized distribution functions for 
quasiparticles. 



systems describing Fermi (Bose) quasiparticles interact- 
ing with spins. The real time formalism is also an alterna- 
tive approach for the analytical continuation method for 
equilibrium problems allowing direct calculations of cor- 
relators whose analytical properties as function of many 
complex arguments can be quite cumbersome. 

A long time ago Keldysh jj8| and Schwinger |59) have 
proposed a novel approach for the description of kinetic 
phenomena in metals. This approach was found especially 
fruitful for normal metals p0|, and, in many recent ap- 

for disordered inter- 



plications, for superconductors 41 
acting (normal or superconducting) electron liquids [[12| 
for example. The previous application of the real-time 
formalism to the quantum theory of Bose-Einstein con- 
densation (BEC) pa] allowed the derivation of a Fokker- 
Planck equation, which describes both kinetic and coher- 
ent stages of BEC. Moreover p(| developed the closed- 
time path integral formalism for aging effects in quantum 
disordered systems being in contact with an environment. 
The Keldysh technique in application to disordered sys- 
tems (see Q-Q and |i7j -Q) has also been recently 
applied to develop a field theory alternative to the previ- 
ously used replica technique. 

To derive the real-time formalism for SU(iV) gener- 
ators we use the path integral representation along the 
closed time Keldysh contour (see Fig j^) . Following the 
standard route |49| we can express the partition function 
of the problem containing SU(iV) generators as a path in- 
tegral over Grassmann variables ipi = (ai t i(j), ai t N(j)) T 
where I = 1,2 stands for upper and lower parts of the 
Keldysh contour, respectively, 

Z/Z = J D^D^exp{iA)/ J D$Dipexp(iAo) (41) 

where the actions A and Ao are taken as an integral along 
the closed-time contour Ct + C T which is shown in Figjg. 
The contour is closed at t = — oo + it since exp(—(3H ) = 

T T cxp Jq Hodr) . We denote the tp fields on upper and 

lower sides of the contour Ct as tpi and ip2 respectively. 
The fields W stand for the contour C T . These fields pro- 
vide matching conditions for ipi t 2 and are excluded from 
final expressions. Taking into account the semi-fermionic 
boundary conditions for generalized Grassmann fields (|2^ ) 
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one gets the matching conditions for i/>i 2 at t = ±00, 

V>l,Jfc(-°o) = ex P ( i7r ~~^~) ^2,o|fc(-0°) 

<alk(+°°)=$UI*(+°°) ( 4 2) 

for fc = 1, |_7V/ 2j and a = 1,...,N. The correlation 
functions can be represented as functional derivatives of 
the generating functional 



Z[rj\ = Zq 1 J DtpDipexp [iA + i j> dt(f]a z tp + ^a z n) 

(43) 

where 77 represents sources and the a z matrix stands for 
"causal" and "anti-causal" orderings along the contour. 

The on-site Green's functions (GF) which are matrices 
of size 2N x 2N with respect to both Keldysh (lower) and 
spin-color (upper) indices are given by 



z[vh, 



r),r}— j-0 



(44) 



To distinguish between imaginary-time (j27]) and real-time 
( f44| ) GF's we use different notations for Green's functions 
in these representations. 

After a standard shift-transformation 0] of fields tp 
the Keldysh GF of free semi-fermions assumes the form 



G%{e) = Gq 



R,a ( 1 - h - ft 



Gi 



A,a ( — ft —ft 



1-fe-fe) ^° VWeW, 

where the retarded and advanced GF's are 

Gi R ' A)a (e) = (e±i6)-\ f t = f^ k \e) (45) 
with equilibrium distribution functions 



/W fc )(e)=T 



1 



_ e e i7r(2*-i)/JV ex p(/3e) + 1 

(46) 

A straightforward calculation of f( N < k > for the case of even 
N leads to the following expression 



N 

£< 

1=1 



(-l) t_1 exp (j3e(N-l)) exp 



M(2k- 1) 
N 



exp(7V/3e) + 1 

(47) 

where k = 1, ...,N/2. 

The equilibrium distribution functions (EDF) f( 2S + 1 - k ) 
for the auxiliary Fermi-fields representing arbitrary 5 for 
SU (2) algebra are given by 



(2S+l,fc)/ £ \ = 



2S+1 



J2 exp {Pe(2S +1-0) exp 



{=1 



in{2k - 1). 
25+1 



exp((25+l)/3e) + (-l) 2 S+ 1 



(48) 



for k = 1, \_S + 1/2J. Particularly simple are the cases 
of 5 = 1/2 and 5 = 1, 



/( 2 ^)(e)=n F (2 e )-z 



1 



2 cosh(/3e) 



(,,1,, 1 , . 3 , . 1— sinh(/3e/2) 

/ <8,1, W = 2^W-2 n ^- iV5 itoh^ (49) 

Here, standard notations for Fermi/Bose distribution func- 
tions n F / B (e) — [exp(/3e) ± 1] _1 are used. 

In general the EDF for half-integer and integer spins 
can be expressed in terms of Fermi and Bose EDF respec- 
tively. We note that since auxiliary Fermi fields introduced 
for the representation of SU(iV) generators do not repre- 
sent the true quasiparticles of the problem, helping only 
to treat properly the constraint condition, the distribution 
functions for these objects in general do not have to be real 
functions. Nevertheless, one can prove that the imaginary 
part of the EDF does not affect the physical correlators 
and can be eliminated by introducing an infinitesimally 
small real part for the chemical potential. In spin prob- 
lems, a uniform/staggered magnetic field usually plays the 
role of such real chemical potential for semi-fermions. 

Let us illustrate the application of the semi-fermionic 
formalism for spin Hamiltonians. As an example we con- 
sider the SU(2) Heisenberg model for 5=1/2 with the 
nearest neighbor interaction 



tut 



<ij> 



(50) 



We start with imaginary-time semi-fermionic descrip- 
tion of the ferromagnetic (FM) state of the Heisenberg 
model (J > 0). We follow the standard procedure devel- 
oped in the original paper of Popov and Fedotov [g5| . Af- 
ter applying the Hubbard-Stratonovich transformation to 
decouple the four-semi-fermion term in (|5(]) by the local 
vector field 3^(t) the effective action is obtained: 

I f drJ2(lFM(<l))- X * q (T)* q (T) (51) 
Jo q 

where i[> T = (V'TV'i) " fields denoting the semi-fermions in 
SU(2) representation of the 5 = 1/2 spin operators, 



hi(q) = Ifm/afm (?) = jj E J ( r y)e iq 



(52) 



and 7m (0) = ZJ > for the FM instability (here Z de- 
notes the number of the nearest neighbors, N stands for 
the number of unit cells). The FM phase transition cor- 
responds to appearance at T = T c of the nonzero average 
(<P(0,0)) which stands for nonzero uniform magnetiza- 
tion, or by another words, corresponds to the Bose con- 
densation of the field <P Z . 

Splitting the field <P Z on the time-independent spatially 
homogeneous (uniform) part and the fluctuating filed <P Z 

$ z (k,u) = A4(/3A0 1/2 (WL,o + <^(k,^) (53) 
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where ^-(q,0) = 2 X g z (q,0) = 2 Xo = 25(5 + 1)0/3 
stands for the on-site spin susceptibility in paramagnetic 
state. 

Let us now consider the Heisenberg model with antifer- 
romagnetic (AFM) sign of the exchange integral (J < 0). 



Fig. 7. First few graphs for the free energy expanded with 
respect to local molecular magnetic field. Solid line stands for 
semi-fermions. Zig-zag line denotes the "condensate" field. 



make it possible to integrate over all semi-fermionic fields. 
As a result, the nonpolynomial effective action can be de- 
rived for the FM Heisenberg model 



A$ f M [1>,*]=Ml>,*]+ 

- f dTY,{lAFM{<l))- l *> q {T)* q {T) 



(59) 



A ef f = A [*] + Tr In (g~ 1 (#* , <S> ± )) 



(54) 



where Q a = — (TV^o-Q, T)tp a (ji 0) stands for the local Green's 
function of semi-fermions. The expansion of the Tr In Q~ l 
with respect to <P fields results in standard Ginzburg- 
Landau functional (see Fig. 7). The effects of molecular 
field are included into zero approximation for GF: 

0°(*w w ) = [iu n + al„M/2]-\ 



In one loop approximation the standard molecular field 
equation can be reproduced 



M = J™(0)tanh(#A/t/2) 



(55) 



The saddle point (mean-field) effective action is given 
by well-known expression 



MM] 



-N 



4/m(0) V V 2 J 



(56) 



and the free energy per spin f (see Fig. 7) is determined 
by standard equation: 



0fo = — In Z s = — hi I 2 cosh ' 



4,1m (0) 



V 2 



(57) 



Calculation of the second variation of A e f / gives rise 
to the following expression 



SA 



eff 



-zE^( k >°) 







$>+(k,c) 



2 cosh 2 (/3 J?) 
tanh(/3f2) 



2f2 



#~(k,w) (58) 



where Q = (gusH + M)/2. For T > T c one easily ob- 
tains the effective static spin-spin interaction equivalent to 
those given by the Random Phase Approximation (RPA) 

T(q, 0) = (*(q, 0)*(-q, 0)> = 2/(q)/(l - 2 X o/(q)), 



and Im(Q) = ZJ < for the AFM instability correspond- 
ing to vector Q = (n, ...,it) (we consider the hypercubic 
lattice for simplicity). In contrast to the FM case, we can 
now represent the longitudinal component of the field <!> z 
as a superposition of the staggered time-independent part 
(" staggered condensate" ) and a fluctuating field 

<P z (k,co) =A^(/3A0 1/2 5 k , Q (L,o + <^(k, W ) (60) 

As a result, the integration over semi-fermionic fields can 
be done explicitely. Introducing two sublattices for ip fields 
one gets 4x4 matrix structure for the semi-fermionic 
Green's functions. Since the AFM instability is associated 
with appearance of a nonzero staggered magnetization Af, 
it is necessary to take into account both "normal" and 
"anomalous" GF determined as follows: 



dre^ T (T T W(k,r)Vv (k,0)> 



lLO n 



(61) 



J*(iw„) = - / dTe^" r (T T Vv(k,T)^(k + Q,0)) = 
Jo 



(62) 



where fl = (Af+gfiBh)/2. Integrating over all semi-fermions 
one obtains the mean-field equation for the staggered mag- 
netization: 



Af = -Iafm(Q) tanh(/W/2) 



<F(k,0)- 



and 



A [Af] = N 



0Af 2 
4,1m (Q) 



In 2 cosh 



0Af 



(63) 



(64) 



After taking into account the second variation of A e f / 
the following expression for the effective action is obtained 
[(see e.g. ||,|f7|): 



^M(k) 







2cosh 2 (/3i7) 
- ^ x (k)**(k, W )*«(k, W )+ 



<P z (k,0)- 
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k,cj 



^M(k) 



2^tanh(/3^) 



4^ 2 



]T<Z>+(k + Q, w ) 



<5 _ (k,w)- 
<2>"(k,cj). (65) 



The application of the Schwinger-Keldysh formalism 
for the Heisenbcrg model is straightforward. Applying the 
semi-fermionic transformation to the partition function 
one obtains the action as an integral along the closed-time 
Keldysh-contour 

A = Aq + A int =A +I dtJ2 J(q)S q (t)S_ q (t) (66) 
J c 



where Aq corresponds to noninteracting semi-fcrmions 





A -/ o *?*( (0S T" 1 (0* i )-'.l-' 



(67) 



We denote ^ q =^X)<i> eiql > ^q^Jq/Jo and apply four - 
component semi-fermionic representation for FM case and 
eight - component representstion with ip T — (ip^ ^Ak+cj) f° r 
the AFM case. Performing the standard Hubbard - Stra- 
tonovich transformation along the Keldysh contour with 
the help of the two-Keldysh-component vector (Bose) field 
one gets 



Ai 



-Tr{$l J- V* q ) + Tr$Q»o-j^) (68) 



Now we integrate out V fields and express the effective 
action in terms of «& fields 

A eff = -irr^J-VSq) + Tr In (G„ 1 + *^) 

where 7 /i = (tr z ± l)/2 acts in Keldysh space. Since in gen- 
eral 3? is a time- and space-dependent fluctuating field 
the partition function (uW cannot be evaluated exactly. 
Nevertheless, when a magnetic instability occurs, we can 
represent the longitudinal component of this field as a su- 
perposition of a uniform (FM) or staggered (AFM) time- 
independent part and a fluctuating field 

<P z ^(q,u)) = condensate + <^(q, w), (69) 

where <£^(q, uo) = (j)^(q_, oj) with the matching conditions 
at t = ±oo 

0f(-oo) = ^(-oo), 0f (+oo) = ^(+oo). (70) 
We expand Tr 1ii(Gq 1 + 4> ll a r ) li ) in accordance with 

°° (_i)n+l 

Tr ln(.. .) = Tr In Go 1 +TrV ^ ^ (Gq^cttT 

n=l 

(71) 

The spectrum of the excitations (FM or AFM magnons) 
can be defined as poles of the transverse GF 



D(x I *)=-t(T c ^+(x,t)0r(O,O)>. 




- 6) 0- 

Fig. 8. Feynman diagrams contributing to dispersion (a) and 
damping (b) of magnons. Solid line denotes semi-fermions. 



The procedure of the calculation of this GF is similar to 
that for a "fcrmionic" GF. Introducing the sources and 
evaluating (fn]) one gets 



R [ 1 + N u N u \ _ A ( N u N u 
l + N u N u ) u ° \ l + N u l + N u 



where the retarded and advanced magnons GF's are 
D R < A {q, uu) = {u>-L>(q)±i6)- 1 , = (exp(/?w) - l)" 1 



The magnon spectrum is determined by the zeros of the 

2 



determinant of J 1 — TI^ (w) (see Fig. 8a) in equilibrium 



w q = J M(1 - !/ q ) \cf, 
for FM magnons and 



(72) 



(73) 



for AFM magnons. The uniform and staggered magnetiza- 
tion are given by equations (^5|) and ( |63| ) correspondingly. 

The magnon damping is defined by four-magnon pro- 
cesses 7T4 1 , shown in Fig8.b. The derivation of the 
kinetic equation and calculation of magnon damping is 
reserved here for a detailed publication. 

We consider now the second possibility to decouple 
the four-fermion term in the Heisenberg model with the 
antiferromagnetic sign of spin-spin interaction which can 
be written in a form equivalent to (pG): 



1 



Hint — — 



2 ^ 



Jiji/jjipji/jUi + - V JijpiPj 



4 ^ 

<ij> 



(74) 



Last term in the Hamiltonian ( f74| ) describes the fluctua- 
tions of semi-fermionic density pi — i/jjipi and therefore 
is irrelevant for our calculations. In contrast to descrip- 
tion of the local correlations achieved with the help of the 
local vector bosonic field we introduce now the bi-local 
scalar bosonic field Aij depending on two sites and re- 
sponsible for inter-site semi-fermionic correlations. Defin- 
ing new coordinates R=(Ri+Rj)/2, r=Ri— Rj and ap- 
plying a Fourier transformation we obtain the effective 
action 



A 



eff 



-Tr(A T Pqi J^a z A Pq2 )+Tr In (Gq 1 - A p7 
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This effective action describes the nonequilibrium quan- 
tum spin-liquid (SL). We confine ourselves to consider 
the uniform phase Q , Jsij] , , Q of Resonant Valence 
Bonds (RVB) in 2D antiferromagnets. It is suitable to 
rewrite the functional in new variables, namely the am- 
plitude A and the phase 0=rA(R), according to formula 

A< ij> (R, r) = A(r) Jj" exp (irA M (R)) (75) 

The exponent in ( f75| ) stands for gauge fluctuations to be 
taken in eikonal approximation. As a result, the effective 
action can be written in continuum limit in terms of the 
gauge fields A a as follows: 

A eff = <j> alt J dkA a (k,t)<rr a/3 Ap(k,t) (76) 

The spectrum of excitations in the uniform SL is defined 
by the zeros of current-current correlation function ]30[| 

^ = Tr(p a pe(G* +(! G« + G^ +q Gf) + S af) f(J p A)) 

in equilibrium |52|,^3| and is purely diffusive (see e.g-fiof ) 

u = iJA\c l \\ 4 = -5>(q)tenh(^) (77) 

We denote G K an off-diagonal element (Keldysh compo- 
nent) of semi-fermionic GF in triangular representation, 
provided that 

G$(e) = -i2Tt5{e ± h)[B 1/2 {(3e) ± isech(/3e)] 

is expressed in terms of a Brillouin function -B1/2 contain- 
ing correct information about occupied states. The equa- 
tion of motion for G K generally constitutes the quantum 
kinetic equation. 

The quantum kinetic equation for nonequilibrium spin 
RVB-liquids can be obtained by taking into account the 
higher order diagrams similarly to Fig. 8b with current-like 
vertices and will be presented elsewhere. 

We discuss now briefly some possible applications of 
the imaginary-time and real-time Schwinger-Keldysh semi- 
fermionic formalism developed for SU(N) Hamiltonians for 
solution of the condensed matter physics problems. The 
Keldysh technique in application to disordered systems at- 
tracts a constant interest (see [^]-|Q) as an alternative 
approach to the replica technique. The main advantage of 
the closed-time contour calculations is an automatic nor- 
malization (disorder independent) of the partition func- 
tion (see |42| ). The application of real-time Schwinger- 
Keldysh approach allows one to study the quantum dy- 
namics of disordered systems being out of equilibrium. We 
note, that the formalism developed in the present paper 
is also a very promising tool for description of a quan- 
tum phase transitions (magnetic, spin-glass etc) in SU(N) 
models (see [^,^9|) Another possible application of the 
semi-fermionic SU(N) representation is the description of 
paramagnet-(ferro) antiferromagnet or paramagnet-spin 
liquid transitions in equilibrium and nonequilibrium stron- 
gly correlated electron systems (see [ |60[|6l| ). The non- 
linear spin waves in strongly correlated local - itinerant 



magnets and the kinetic properties of the nonequilibrium 
spin liquid are also possible problems to be considered 
with the method proposed. The third interesting exam- 
ple of the application of the semi-fermionic formalism is 
the Kondo systems |32|, for example the Kondo lattice 
model usually used for interpretation of an exotic prop- 
erties of heavy-fermion compounds or the nonequilibrium 
Kondo-systems in semiconducting hetero-structures (see 
e.g. J6^]-@|). The main advantage of the semi-fermionic 
representation in applications to the strongly correlated 
systems in comparison with another methods is that the 
local constraint is taken into account exactly and the usual 
Feynman diagrammatic codex for the composite itinerant- 
local compound is applicable. 

Summarizing, we constructed a general scheme for the 
semi-fermionic representation for generators of the SU(iV) 
algebra. A representation for the partition function is found 
both in imaginary and real time. The approach developed 
leads to the standard diagram technique for Fermi opera- 
tors, although the constraint is taken into account rigor- 
ously. The method proposed allows to treat SU(7V) gen- 
erators on the same footing as Fermi and Bose systems. 
The technique derived can be helpful for the description 
of quantum systems in the vicinity of a quantum phase 
transition point and for nonequilibrium systems. 

This work was started on the workshop " Strongly cor- 
related systems" in the ICTP in Trieste. One of the au- 
thors (MNK) is grateful to many participants of this work- 
shop and especially N.Andrei, A.Tsvelik and A.Protogenov 
for fruitful discussions. We thank also D.Aristov and A.Lu- 
ther for valuable comments and interest to our work. The 
work is supported by the DFG under SFB 410, by the 
Evangelisches Studienwerk Villigst (HF) and by an Alexan- 
der von Humboldt fellowship (MNK). 
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